Estimating hitting probabilities of an interacting particle system on a graph
نویسنده
چکیده
We consider the problem of estimating hitting probabilities related to a class of interacting particle systems. These systems, in which two types of particles — ‘electrons’ and ‘holes’ — move on a graph, are simplified versions of models describing charge transport in disordered materials. The probability of interest is the probability that an electron reaches a certain region of the graph before colliding with a hole. We provide a detailed description of our model and explain how it can be simulated. Next, we give a brief introduction to importance sampling, which we use to improve the efficiency of our estimators. To our knowledge, importance sampling has not yet been used to estimate probabilities related to interacting particle systems. We describe how importance sampling can be used to improve the efficiency of estimators of hitting time probabilities involving discrete time Markov chains. We then use importance sampling to estimate the probability that we are interested in. In doing so, we observe that there are a number of complexities that arise when working with interacting particle systems. We describe some simple heuristics for implementing importance sampling. These heuristics make minimal changes to the probability measure under which the original system evolves. We consider a specific example of our problem and investigate the effectiveness of the importance sampling approach. We show that our estimators outperform standard Monte Carlo estimators. Finally, we describe possible future work, which includes a more sophisticated importance sampling approach that uses ‘locally optimal’ changes of measure.
منابع مشابه
A functional central limit theorem for interacting particle systems on transitive graphs
A nite range interacting particle system on a transitive graph is considered. Assuming that the dynamics and the initial measure are invariant, the normalized empirical distribution process converges in distribution to a centered di usion process. As an application, a central limit theorem for certain hitting times, interpreted as failure times of a coherent system in reliability, is derived.
متن کاملGenetic genealogical models in rare event analysis
We present in this article a genetic type interacting particle systems algorithm and a genealogical model for estimating a class of rare events arising in physics and network analysis. We represent the distribution of a Markov process hitting a rare target in terms of a Feynman–Kac model in path space. We show how these branching particle models described in previous works can be used to estima...
متن کاملSequential Monte Carlo simulation for the estimation of small reachability probabilities for stochastic hybrid systems
The problem of estimating the probability that a system reaches a given set within some time horizon is considered. Standard Monte Carlo methods for reachability probability estimation do not require specific assumptions on the system under consideration. However, they are computationally demanding when the probability to be estimated is small. An Interacting Particle System (IPS) approach has ...
متن کاملRare Event Simulation Using Reversible Shaking Transformations
We introduce random transformations called reversible shaking transformations which we use to design two schemes for estimating rare event probability. One is based on interacting particle systems (IPS) and the other on time-average on a single path (POP) using ergodic theorem. We discuss their convergence rates and provide numerical experiments including continuous stochastic processes and jum...
متن کاملMean-Field Limit of a Stochastic Particle System Smoothly Interacting Through Threshold Hitting-Times and Applications to Neural Networks with Dendritic Component
In this article we study the convergence of a stochastic particle system that interacts through threshold hitting times towards a novel equation of McKeanVlasov type. The particle system is motivated by an original model for the behavior of a network of neurons, in which a classical noisy integrate-and-fire model is coupled with a cable equation to describe the dendritic structure of each neuron.
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 2015